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Abstract 

In this article, we use a mechanism first introduced by Herman, 
Marco, and Sauzin to show that if a Gevrey or analytic perturba- 
tion of a quasi-convex integrable Hamiltonian system is not too small 
with respect to the number of degrees of freedom, then the classical 
exponential stability estimates do not hold. Indeed, we construct an 
unstable solution whose drifting time is polynomial with respect to the 
inverse of the size of the perturbation. A different example was already 
given by Bourgain and Kaloshin, with a linear time of drift but with a 
perturbation which is larger than ours. As a consequence, we obtain a 
better upper bound on the threshold of validity of exponential stability 
estimates. 



1 Introduction 

1. Consider a near- integrable Hamiltonian system of the form 

--hii) + f{e,i) 




with angle-action coordinates (0,1) G T" x R", and where / is a small 
perturbation, of size e, in some suitable topology defined by a norm [ . |. 

If the system is analytic and h satisfies a generic condition, it is a remark- 
able result due to Nekhoroshev ( |Nek77j . |Nek79| ) that the action variables 
are stable for an exponentially long interval of time with respect to the 
inverse of the size of the perturbation: one has 

\Iit) -Io\< cie\ \t\ < C2 exp(c3e~"), 

for some positive constants ci,C2,C3,a, 6 and provided that the size of the 
perturbation e is smaller than a threshold Eq. Of course, all these constants 
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strongly depend on the number of degrees of freedom n, and when the latter 
goes to infinity, the threshold Eq and the exponent of stability a go to zero. 

More precisely, in the case where h is quasi-convex and the system is 
analytic or even Gevrey, then we know that the exponent a is of order 
n"^ and this is essentially optimal (see |LN92j . |Pos93j . [MSn2j . |LMn5j 
[Zha09) and [BMlOj for more information on the optimality of the stability 
exponent). 

2. This fact was used by Bourgain and Kaloshin in |BK05] to show that if 
the size of the perturbation is 



then there is no exponential stability: they constructed unstable solutions 
for which the time of drift is linear with respect to the inverse of the size of 
the perturbation, that is 

More precisely, in the first part of [BK05j . Bourgain proved this result for 
a specific example of Gevrey non-analytic perturbation of a quasi-convex 
system, then for an analytic perturbation he obtained a time t„ ~ ^n^~^-, 
for any c > 0. In the second part of |BK05] . using much more elaborated 
techniques (especially Mather theory), Kaloshin proved the above result for 
both Gevrey and analytic perturbation and for a wider class of integrable 
Hamiltonians, including convex and quasi-convex systems. 

Their motivation was the implementation of stability estimates in the 
context of Hamiltonian partial differential equations, which requires to un- 
derstand the relative dependence between the size of the perturbation and 
the number of degrees of freedom. Their result indicates that for infinite 
dimensional Hamiltonian systems, Nekhoroshev's mechanism does not sur- 
vive and that "fast diffusion" should prevail. Of course, in their example, 
one cannot simply take n = oo as the size of the perturbation £n ~ e~" goes 
to zero and the time of instability r„ ~ goes to infinity exponentially fast 
with respect to n. A more precise interpretation concerns the threshold of 
validity Eq in Nekhoroshev's theorem: it has to satisfy 

eo « e"", 

and so it deteriorates faster than exponentially with respect to n. 

3. As was noticed by the authors in [BK05j . their examples share some 
similarities with the mechanism introduced by Herman, Marco and Sauzin 
in [MS02| (see also [LM05] for the analytic case). In this present article, we 
use the approach of [MS02j and [LM05j to show, using simpler arguments 
than those contained in [BK05j . that if the size of the perturbation is 

e„ ~ g-nln(nlnn)^ 
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then it is still too large to have exponential stability: we will show that one 
can find an unstable solution where the time of drift is polynomial, more 
precisely 

\Iirn) - Io\ ~ 1, T„ ~ e-"". 

As in the first part of |BK05j . we will construct specific examples of Gevrey 
and analytic perturbations of a quasi-convex system. We refer to Theo- 
rem [27T] and Theorem 12.31 below for precise statements. Hence one can infer 
that the threshold £q in Nekhoroshev's theorem further satisfies 

eo « e-" '"("i'^") « e"", 

and this gives another evidence that the finite dimensional mechanism of 
stability cannot extend so easily to infinite dimensional systems. Let us 
point out that our time of drift is worst than the one obtained in (BK05], 
but this stems from the fact that the size of our perturbation is smaller than 
theirs and so it is natural for the time of instability to be larger. Moreover, 
our exponent n in the time of drift can be a bit misleading since in any cases, 
that is even for a linear time of drift, t„ goes to infinity exponentially fast 
with n, so such results do not apply at all to infinite dimensional Hamiltonian 
systems. A natural question, which was raised by Marco, is the following. 

Question 1. Given e > arbitrarily small, construct an £fi-perturbation of 
an integrable system having an unstable orbit with a time of instability Tn 
such that 

lim En = e, lim r„ < +oo. 

n— >-+oo n— >-+oo 

We believe that one can give a positive answer to this question, by using 
a more clever construction. However, even if one can formally let n goes 
to infinity, by no means this implies the existence of an unstable solution 
for an infinite-dimensional Hamiltonian systems, which is a very difficult 
problem (see |CKG"'"10 and [GGIOJ for related results in some examples of 



Hamiltonian partial differential equations). 

2 Main results 
2.1 The Gevrey case 

1. Let us first state our result in the Gevrey case. Let n > 3 be the number 
of degrees of freedom, and given i? > 0, let i? = be the open ball of 
around the origin, of radius R > with respect to the supremum norm | . |, 
and B its closure. 

The phase space is T" x i?, and we consider a Hamiltonian system of the 
form 

H{e,i) = h{i) + f{e,i), {e,i)eT^xB. 
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Our quasi-convex integrable Hamiltonian h is the simplest one, namely 

h{I) = + . . . + + / = (/^, ...,/„)€ i?. 

Let us recall that, given a > 1 and L > 0, a function / S (^oo^ip" x i?) is 
(q, L)-Gevrey if, using the standard multi-index notation, 

|/kL= Yl i'''"(^!)-"|5'/lc0(T"xiJ)<O°- 

The space of such functions, with the above norm, is a Banach space that 
we denote by G'^-^(T" x B). One can see that analytic functions correspond 
exactly to a = 1. 

2. Now we can state our theorem. 

Theorem 2.1. Let n>3, R>1, a>l and L > 0. Then there exist posi- 
tive constants c, 7, C and hq G N* depending only on R, a and L such that for 
any n > uq, the following holds: there exists a function fn € C^'^iT^ x B) 
with En = \fn\a,L Satisfying 

„-2(n-2) ln(4nln2n) ^ ^ ^ „ „-2(n-2) ln(n In 2n) 

such that the Hamiltonian system Hn = h + fn has an orbit {9{t),I{t)) for 
which the estimates 

|/(r„)-/o| > 1, rn<c(— 

hold true. 

As we have already explained, this statement gives an upper bound on 
the threshold of applicability of Nekhoroshev's estimates, which is an im- 
portant issue when trying to use abstract stability results for "realistic" 
problems, for instance for the so-called planetary problem (see [Nie96] ) . 

So let us consider the set of Gevrey quasi-convex integrable Hamiltonians 
n = n{n,R,a,L,M,m) defined as follows: h e H if h G G^'^iB) and 
satisfies both 

\d''h{I)\<M, 1 < |A:i| + ••• + |A:„| < 3, 

and 

V/ £ B,\/ve M", Vh{I).v = ^ S/^h{I)v.v > m\v\^. 

From Nekhoroshev's theorem (see [ MS02j for a statement in Gevrey classes), 
we know that there exists a positive constant eo(^) = £Q{n,R,a,L,M,m) 
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such that the fohowing holds: for any h £ 7i, there exist positive constants 
ci , C2 , C3 , a and b such that if 

/gG"'^(T"x;b), |/U,L<eo(^), 

then any solution {9{t),I{t)) of the system H = h + f , with /(O) S Bji/2, 
satisfies 

- /ol < Cle^ |t| < C2 exp(c3e-"). 
Then we can state the following corollary of our Theorem 12.11 

Corollary 2.2. With the previous notations, one has the upper bound 

eo(^)<e-2("-2)ln(4„ln2n)_ 

This improves the upper bound £o{'H) < e~" obtained in |BK05j for 
Gevrey functions. 

2.2 The analytic case 

3. Let us now state our result in the analytic case. Here B = Bn is still 
the open ball of M" around the origin, of radius -R > with respect to the 
supremum norm, and we will also consider a Hamiltonian system of the form 

H{e,i) = h{i) + f{ej), (0,/)eT"xs, 

where 

m = + ■■■+ Il-l) +In, 1= ih, ...,/„) G 

Given p > 0, let us introduce the space Ap{T^ x B) of bounded real-analytic 
functions on T" x S admitting a bounded holomorphic extension to the 
complex neighbourhood 

Vp = y,(T" xB) = {{0,1) G (C7Z") X I \I{e)\ < p, d{I,B) < p}, 

where I{9) is the imaginary part of 6 and the distance d is associated to the 
supremum norm on C". Such a space Ap{T^ x B) is obviously a Banach 
space with the norm 

l/lp = l/lco(y,) = sup \fiz)\, f G ^,(T" X B). 

zeVp 

Furthermore, for bounded real-analytic vector-valued functions defined on 
X B admitting a bounded holomorphic extension to Vp, we shall extend 
this norm componentwise (in particular, this applies to Hamiltonian vector 
fields and their time-one maps). 

4. Now we can state our theorem. 
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Theorem 2.3. Let n > 4, i? > 1, and cr > 0. Then there exist positive 
constants p, 7, C and uq £ N* depending only on R and a, and a constant 
Cn that may also depends on n, such that for any n > no, the following holds: 
there exists a function fn G ^p(T" x B) with En = \ fn\p satisfying 



such that the Hamiltonian system Hn = h + fn has an orbit (9{t),I{t)) for 
which the estimates 



In the above statement, the constant a has to be chosen sufficiently small 
but independently of the choice of n and R (see Proposition 14. ip . Moreover, 
the theorem is slightly different than the one in the Gevrey case, since there 
is a constant c„ depending also on n: this comes from the use of suspension 
arguments due to Kuksin and Poschel ( |Kuk93j . |KP94j ) in the analytic case, 
which are more difficult than in the Gevrey case. 

Here we can also define a threshold of validity in the Nekhoroshev the- 
orem eo(T-L) = eo{n, R, p, M,m) and state the following corollary of our 
Theorem 12.31 

Corollary 2.4. With the previous notations, one has the upper bound 



This improves the upper bound eo(^) < obtained in |BK05] for 
analytic functions. For concrete Hamiltonians like in the planetary problem, 
the actual distance to the integrable system is essentially of order 10~^, 
hence the above corollary yields the impossibility to apply Nekhoroshev's 
estimates for n > 3. 

2.3 Some remarks 

5. Theorem 12.11 and Theorem 12.31 are obtained from the constructions in 
|MS02j and |LM05j . but one has to choose properly the dependence with 
respect to n of the various parameters involved. 

As the reader will see, we will use only rough estimates leading to the 
factor n in the time of instability: this can be easily improved but we do 
not know if it is possible in our case (that is with a perturbation of size 
g-nin(ninn)^ to obtain a linear time of drift. 

Let us note also we have restricted the perturbation to a compact subset 
of A" = T" X just in order to evaluate Gevrey or analytic norms. In fact, 



-2(n-3)ln(4nln2n) ^ ^ -2{n-3) ln(ra In 2ra) 




hold true. 



eo{n) < e 



2(n-3) ln(4n In 2n) 
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in both theorems the Hamiltonian vector field generated by Hn = /i + is 
complete and the unstable solution {6{t),I{t)) satisfies 

lim \I{t) — Iq\ = +00, 

t— i-±oo 

which means that it is bi-asymptotic to infinity. 

6. li is important to note that our approach leads, as in the first part of 
|BK05] ■ to results for an autonomous perturbation of a quasi-convex inte- 
grable system (or, equivalently, for a time-dependent time-periodic pertur- 
bation of a convex integrable system). In the second part of |BK05j . for 
a class of convex integrable systems, Kaloshin was able to reduce the case 
of an autonomous perturbation to the case of a time-dependent time peri- 
odic perturbation, partly because of his more general (but more involved) 
approach. We could have tried to apply his general arguments to our case, 
but for simplicity we decided not to pursue this further. 

7. In this text, we will have to deal with time-one maps associated to 
Hamiltonian fiows. So given a function H, we will denote by the time-t 
map of its Hamiltonian flow and by = the time-one map. We shall 
use the same notation for time-dependent functions H, that is will be 
the time-one map of the Hamiltonian isotopy (the flow between t = and 
t = I) generated by H. 

3 Proof of Theorem 12.1 

The proof of Theorem 12.11 is contained in section 13.21 but flrst in section 13.11 
we recall the mechanism of instability presented in the paper |MS02j (see 
also (MSOlj). 

This mechanism has two main features. The flrst one is that it deals with 
perturbation of integrable maps rather than perturbation of integrable flows, 
and then the latter is recovered by a suspension process. This point of view, 
which is only of technical matter, was already used for example in [Dou88] 
and offers more flexibility in the construction. The second feature, which is 
the most important one, is that instead of trying to detect instability in a 
map close to integrable by means of the usual splitting estimates, we will 
start with a map having already "unstable" orbits and try to embed it in 
a near-integrable map. This will be realized through a "coupling lemma", 
which is really the heart of the mechanism. 

As we will see, the construction offers an easy and very efficient way of 
computing the drifting time of unstable solutions, therefore avoiding all the 
technicalities that are usually required for such a task. 
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Figure 1: Drifting point for the map i/jg 



3.1 The mechanism 

1. Given a potential function [/ : T — )• M, we consider the following family 
of maps ipg : A ^ A defined by 

il;g{e,I) = {0 + qI,I-q-'U'{e + qI)), (9,1) e A, (1) 

for g E N*. If we require U'{0) = —1, for example if we choose 

U{e) = -(27r)-^ sin(27r0), U'{9) = - cos(2^0), 

then it is easy to see that 'ipg{0,0) = (0, g^^) and by induction 

^^;'^{0,0) = {0,kq-') (2) 

for any A; E Z (see figured]). After q iterations, the point (0, 0) drifts from the 
circle / = to the circle 1 = 1 and it is bi-asymptotic to infinity, in the sense 
that the sequence {ipq{0, 0))^^^ ^® ^^^^ contained in any semi-infinite annulus 
of A. Clearly these maps are exact-symplectic, but obviously they have no 
invariant circles and so they cannot be "close to integrable". However, we 
will use the fact that they can be written as a composition of time-one maps, 

V;, = ^l-'U o ($1^' o . . . o $5^') = ^l-'U o ($1^')' , (3) 

to embed ijjq in the g'^^-iterate of a near-integrable map of A", for n > 2. To 
do so, we will use the following "coupling lemma", which is easy but very 
clever. 
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Lemma 3.1 (Herman-Marco-Sauzin). Let m,m' > 1, F : A™ A*" and 

G : A'"' A*"' two maps, and f : A™ M and g : A™' -> R tu^o Hamil- 
tonian functions generating complete vector fields. Suppose there is a point 
a G A"^ which is q-periodic for G and such that the following "synchronisa- 
tion" conditions hold: 

5(a) = 1, dg{a)=0, r7(G^(a)) = 0, dgiG''{a)) = 0, (S) 

for 1 < k < q — 1. Then the mapping 

q, ^ o{FxG): A™+'"' — > A'"+™' 

is well-defined and for all x € A*", 

■^'i{x,a) = oF'?(x),a). 

The product of functions acting on separate variables was denoted by (8", 

1. e. 

f(g)g{x,x') = f{x)g{x') x e A™, x' G A™'. 
Let us give an elementary proof of this lemma since it is a crucial ingredient. 

Proof. First note that since the Hamiltonian vector fields Xf and Xg are 
complete, an easy calculation shows that for all x G A*", x' € A™ and t G R, 
one has 

$f®^(x,x') = (l>f("')^(x),${(")^(x')) (4) 

and therefore Xf^g is also complete. Using the above formula and condi- 
tion dH]), the points {F^{x), G^{a)), for 1 < k < q-1, are fixed by ^^^^ and 
hence 

^''-\x,a) = {F'i-\x),G'i-\a)). 
Since a is q-periodic for G this gives 

*'?(x,a) =$-^®^'(F<?(x),a), 

and we end up with 

using once again ^ and Q. □ 

Therefore, if we set m = 1, F = $2^1 and / = q^^U in the coupling 
lemma, the g^'^-iterate ^''^ will leave the submanifold A x {a} invariant, 
and its restriction to this annulus will coincide with our "unstable map" 
Hence, after q'^ iterations of the /i-component of the point ((0, 0), a) G A^ 
will move from to 1. 

2. The difficult part is then to find what kind of dynamics we can put 
on the second factor to apply this coupling lemma. In order to have a 
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continuous system with n degrees of freedom at the end, we may aheady 
choose m' = n — 2 so the couphng lemma will give us a discrete system with 
n — 1 degrees of freedom. 

First, a natural attempt would be to try 

1 r2 I I 1 r2 
G=Gn = «>2-f2+-+2^n-l. 

Indeed, in this case 

FxGn = $i^2+-+i^n-i = 

where 

and the unstable map ^ given by the coupling lemma appears as a pertur- 
bation of the form ^' = ^"o^'^, with u = f (^g. However, this cannot work. 
Indeed, for j G {2, . . . ,n — 1}, one can choose a pj-periodic point a'--'-* € A 

i/2 

for the map $ 2 i , and then setting 

a„ = (a(2), . . . , a("-i)) G A"-^, q^=p^. ..p^_,, 

the point a„ is g„-periodic for Gn provided that the numbers pj are mutually 
prime. One can easily see that the latter condition will force the product 
Qn to converge to infinity when n goes to infinity. So necessarily the point 
a„ gets arbitrarily close to its first iterate Gn{cin) when n (and therefore g„) 
is large: this is because g^-periodic points for Gn are equi-distributed on 
g„-periodic tori. As a consequence, a function gn with the property 

gn{an) = 1, 5n(G„(a„)) = 0, 

will necessarily have very large derivatives at a„ if qn is large. Then as the 
size of the perturbation is essentially given by 

\f ®9n\ = \qn^U ®gn\ = \qn^\\gnV 

one can check that it is impossible to make this quantity converge to zero 
when the number of degrees of freedom n goes to infinity. 

3. As in [MSn2] . the idea to overcome this problem is the following one. We 
introduce a new sequence of "large" parameters Nn G N* and in the second 
factor we consider a family of suitably rescaled penduli on A given by 

Pn{02j2) = \ll + N-^V{e2), 

where V{6) = — cos27r0. The other factors remain unchanged, so 

Gn = $i(^l+^l + -+-fn-l) + ^n'V'{02)_ 
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In this case, the map ^ given by the couphng lemma is also a perturbation 
of but of the form ^ = ^^o xJ/'^+^j with v = N~'^V. But now for the map 
Gn, due to the presence of the pendulum factor, it is now possible to find 
a periodic orbit with an irregular distribution: more precisely, a g^-periodic 
point a„ such that its distance to the rest of its orbit is of order A'^^, no 
matter how large g„ is. 

4. Let us denote by {pj)j>o the ordered sequence of prime numbers and let 
us choose Nn as the product of the n — 2 prime numbers {pn+a, • • • ■,P2n}-, 
that is 

Nn = Pn+3Pn+i ' ' ' P2n & ^* ■ 

Our goal is to prove the following proposition. 

Proposition 3.2. Let n > 3, a > 1 and Li > 0. Then there exist a 
function Qn G G"'^^ (A"'"^), a point a„, G A"~^ and positive constants c\ and 
C2 depending only on a and Li such that if 



Cl 



^ C2(n-2)p- 



then an is qn-periodic for Gn and (gn, Gn, an,qn) satisfy the synchronization 
conditions 

gn{an) = 1, dgn{an) = 0, gn{G^{an)) = 0, dgn{G^{an)) = 0, 
for 1 < k < qn — 1- Moreover, the estimate 

holds true. 

The rest of this section is devoted to the proof of the above proposition. 
Note that together with the coupling lemma (Lemma 13. ip . this proposition 
easily gives a result of instability (analogous to Proposition 2.1 in [MS02j) 
for a discrete system which is a perturbation of the map ^>'', but we prefer 
not to state such a result in order to focus on the continuous case. 

5. We first consider the simple pendulum 

P{9j) = h^ + V{e), (6,1) e A. 

With our convention, the stable equilibrium is at (0, 0) and the unstable one 
is at (0, 1/2). Given any M £ N*, there is a unique point b^^ = {OJm) which 
is M-periodic for (this is just the intersection between the vertical line 
{0} X R and the closed orbit for the pendulum of period M). One can check 
that Im £]2,3[ and as M goes to infinity, (0,/m) tends to the point (0,2) 
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—a cr 

Figure 2: The point and its iterates 

which belongs to the upper separatrix. Since Pn{9,I) = + N~'^V{6), 
then one can see that 

where SniO,I) = {9,NnI) is the rescahng by Nn in the action components. 
Therefore the point b^^ = (0, N~^Im) is gn-periodic for for g„ = NnM. 
Let ($f be the flow of the pendulum, and 

$f(o,/M) = {eM{t),iM{t)). 

The function 9M{t) is analytic. The crucial observation is the following 
simple property of the pendulum (see Lemma 2.2 in |MS02j for a proof). 

Lemma 3.3. Let a = -\ + ^ arctane'' < \. For any M e N*, 

Ouit) i [-cj,o-], 

/ort G [1/2,M- 1/2]. 

Hence no matter how large M is, most of the points of the orbit of 
b^ € A will be outside the set {-a < 6 < a} x R (see figured]). The 
construction of a function that vanishes, as well as its first derivative, at 
these points, will be easily arranged by means of a function, depending only 
on the angle variables, with support in {—a < 9 < a}. 

As for the other points, it is convenient to introduce the function 

TM : [-a, a] ^] -l/2,l/2[ 
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which is the analytic inverse of 9m- One can give an expUcit formula for 
this map: 



TM 



(0)= [ 
Jo 



l\, — 4 sin^ Tup 



In particular, it is analytic and therefore it belongs to G"''^^ ([— cj, o"]) for 
a > 1 and Li > 0, and one can obtain the following estimate (see Lemma 
2.3 in [MSn2j for a proof). 

Lemma 3.4. For a > 1 and Li > 0, 

A = sup |rA/|Q,Li < +00. 

AfeN* 

Note that A depends only on a and Li. Under the action of tm, the 
points of the orbit of whose projection onto T belongs to {—a < ^ < o"} 
get equi-distributed, and we can use the following elementary lemma. 

Lemma 3.5. For p G N*, the analytic function 7]p iT ^M. defined by 

r]p{e) = -^cos27r/0 
\^ 1=0 J 

satisfies 

r?p(0) = 1, ??;(0) = 0, vp{k/p) = v'p{k/p) = 0, 
for i < k < p — 1, and 

1 77 I r < p2aLi{27rp)Q 
\'lp\ci,Li '-■ 

The proof is trivial (see |MS02j . Lemma 2.4). 
6. We can now pass to the proof of Proposition 13.21 

Proof of Proposition \3.2l For a > 1 and Li > 0, consider the bump function 
^a,Li £ G"'^i(T) given by Lemma fA. II (see Appendix . 

We choose our function gn G G"'^^(A"^^), depending only on the angle 
variables, of the form 



where 



and 

5i*^(^0= 3<^<n-L 

Let us write 



Cl 



03 _i_ 
Q!,(4ct) a Li 



a,(4o-) a Li 
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Now we choose our point a„ = {an \ • • • , fln" ^- We set 

and 

a« = (0,p;|i+,), 3<i<n-l. 
Let us prove that (Xn is (/^-periodic for Gn^ ^Vg can write 

Since Pn+i, ■ ■ ■ ,P2n are mutually prime, the point {an \ . . . ,a"n G A"~^ 
is periodic for G, with period 

N'n = Pn+i ■■■P2n- 

By construction, the point = 6*^" € A is periodic for with period 
Qn = NnMn, where 

Nn = Pn+3Pn+A ' ' ' P2n- 

This means that a„ is periodic for the product map G„, and the exact period 
is given by the least common multiple of g„ and N'^. S ince divides (i^i^ 
the period of an is qn- 

Now let us show that the synchronization conditions dS]) hold true, that 

is 

Qnian) = 1, dgn{an) = 0, gn{Gl{an)) = 0, dgn{Gi{an)) = 0, 
for 1 < /c < g„ — 1. Since (/?a,Li(0) = 1, then 

9n(a„) = <7l'H0)...5i"-^)(0) = l 
and as ip'^ l^(0) = 0, then 

dgn{an) = 0. 

To prove the other conditions, let us write G^(a„) = {9k,Ik) G A""^, for 
l<k<qn-l. 



^2) (2) 

If 0^ does not belong to ] — cr, cj[, then gn and its first derivative vanish 



(2) 

at Oi because it is the case for w , , i , so 

^ a,{4:a) aL 



9n{6k) = dgn{Ok) = 0. 

(2) 

Otherwise, if —a < Oj^ < a, one can easily check that 

< k < 



and therefore 

3(2)^ = A 
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while 

= 3 < i < n - 1. 

Pn+i+l 

If A'^^ = Pn+i • • • P2n divides k, that is k = k'N^ for some k' G Z, then 

= 1^ = 

. (2) 

and therefore, by Lemma \'6.5\ Vpn+s vanishes with its differential at , and 
SO does gn ■ Otherwise, does not divide k and then, for 3 < i < n — 1, 

(2) 

at least one of the functions r/p^_|^^^. vanishes with its differential at , and 
so does gn^ . Hence in any case 

gniOk) = dgn{0k) = 0, l<k<qn-l, 

and the synchronization conditions ([S]) are satisfied. 

Now it remains to estimate the norm of the function g^- First, using 
Lemma lA.21 one finds 

|5n|a,Li ^ 

which by definition of ci gives 

n + 3 ° '^M-n \a,Li |^p„+4 I^P2n \a,Ll ■ 

Then, by definition of A (Lemma l3.4p and using Lemma fA.31 (with Ai = A^ 
SO using Lemma [331 and setting C2 = 2a sup |a«,Li| (27r)Q, this gives 



09 _1 
a,(4(T) a Li 



1 \Vpn+3 ° '''M„\a,Li\'np„+4\a,Li 
a, (4a) a Li 



I I <; 2Q:(A- + (n-3)Li)(27rp2n)^ 

\yn\a,Li <-lc 

^ ^ ^2asup|AQ ,Li|(n— 2)(27rp2n)^ 

Finally, by definition of M„ we obtain 
and as q„, = NnMn, we end up with 

Qn'^\9n\a,Li < . 

This concludes the proof. □ 
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3.2 Proof of Theorem 



7. In the previous section, we were concerned with a perturbation of the 
integrable diffeomorphism which can be written as o So now 

we will briefly describe a suspension argument to go from this discrete case 
to a continuous case (we refer once again to [MS02] for the details). 

Here we will make use of bump functions, however the process is still 
valid, though more difficult, in the analytic category, (see for example |Dou88j 
or [KP94| ) . The basic idea is to find a time-dependent Hamiltonian function 
on A" such that the time-one map of its isotopy is $"0$'*+^, or, equivalently, 
an autonomous Hamiltonian function on A""*"^ such that its first return map 
to some 2n-dimensional Poincare section coincides with our map o $^+^. 

Given a > 1 and L > 1, let us define the function 

4>Q,L ^ (^j ^^,1^ ^a,L, 

where ipa,L is the bump function given by Lemma lA. II If 4>o{t) = 4>a,L i) 
and <^i(t) = 4>a,L{t — I), the time-dependent Hamiltonian 

H*{e, I, t) = {h{i) + v{e)) ® Mt) + ® Mt) 

clearly satisfies 

But as u and v go to zero, H* converges to h ® (pQ rather than h. How- 
ever, using classical generating functions, it is not difficult to modify the 
Hamiltonian in order to prove the following proposition (see Lemma 2.5 in 
[MSQ2j). 

Proposition 3.6 (Marco-Sauzin). Let n > 1, i? > 1, a > 1, Li > and 

L > satisfying 

= L-(l + (L- + ii+l/2)|</.,,iU,i). (6) 

Ifun,Vn G G"''^i(T"-i), there exists U G G°'-^(T" x^B), independent of the 
variable In, such that if 

Hn{e, I) = ^(l2 + . . . + +In + fn{0, 1), {9, I) G A", 

for any energy e G M, the Poincare map induced by the Hamiltonian flow of 
Hn on the section {9n = 0} n H~^{e) coincides with the diffeomorphism 

o <I)'^+^". 

Moreover, one has 

SUp{|u„|cO, |fn|co} < \ fn\a,L < C3 SUp{ |q:,Li , l^n U,Li }, (7) 

where C3 = 2\(j)a,L\a,L depends only on a and L. 
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8. Now we can finally prove our theorem. 

Proof of Theorem \2.1\ Let ii > 1, a > 1 and L > 0, and choose Li satisfying 
the relation The constants ci and C2 of Proposition 13.21 depend only on 
a and Li, hence they depend only on i2, a and L. 
We can define u„,t;„ G c^.-C'i (T"-i) by 

where U{9i) = — (27r)^^ sin 27r0i, ^(6*2) = — cos27r02 (so Vn is formally 
defined on T but we identify it with a function on T""""*^) and is the 
function given by Proposition 13. 2[ Let us apply Proposition 13.61 there exists 
fn G (7°'-^ (IT" X i?), independent of the variable /„, such that if 

Hn{e,I) = ^{lf + ---+ll_,)+In + fn{0,I), {0,1) €A^, 

for any energy e € M, the Poincare map induced by the Hamiltonian flow of 
H on the section {6n = 0} fl H~^{e) coincides with the diffeomorphism 

Let us show that our system Hn has a drifting orbit. First consider its 
Poincare section defined by 

with 

f^ = q-'U, F = c^H, Gn = ^WI+^i+-+^'n-D+N^'vm_ 
By Proposition 13.21 we can apply the coupling lemma (Lemma 13. ip , so 

^'^((0,0),a„) = (c^J:oF^"(0,0),a„). 
Then, using ([3]), observe that 

$/" o = ^''"''^ o (^^i^^Y" = iPg^, 

so 

^^"((0,0), an) = (($^"oF'?")''"(0,0),a„) 
= ((0,1), 

where the last equality follows from 1^. Hence, after iterations, the Ii- 
component of the point Xn = ((0,0), a„) € A"~^ drifts from to 1. Then, 
for the continuous system, the initial condition (x„,t = 0, I„ = 0) in A" 
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gives rise to a solution {x{t),t, In{t)) = {x{t) , 0n{t) , In{t)) of the Hamiltonian 
vector field generated by Hn such that 

x{k) = kez. 

So after a time r„ = q'^, the point (x„,(0,0)) drifts from to 1 in the 
/i-direction, and this gives our drifting orbit. 

Now let £n = \fn\a,Li be the size of our perturbation. Using the esti- 
mate (l5|) and ([7]) one finds 



(8) 



By the prime number theorem, pn is equivalent to nlnn, so there exists 
no € N* such that for n > uq, one can ensure that 

P2n/4 < Pn+i <P2n, S < i < H, 

which gives 

(W4r"'<A^n<P^„"', Nr' < P2n < ^Nr' . (9) 

We can also assume by the prime number theorem that for n > no, one has 
2n In 2n < p2n < 2(2n In 2n) = 4n In 2n. (10) 
From the above estimates Q and pop one easily obtains 

^(n-2)\n{2-^n\n2n) ^ jy ^ g(?i-2) ln(4n In 2n) ^-^-^-^ 

and, together with ([8]), one finds 

g-2(n-2) ln(4nln2n) ^ g ^ ^^^-2{n-2)ln{2~^nln2n) ^-^2^ 

Finally it remains to estimate the time r„. First recall that 



M„ = 2 



and with 



Hence 



a(n — 2) 



Then using (|lip we have 



Nn^''-'^ < (4nln2n); 
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and from ([8|) we know that 

^ if-" ' 

SO we obtain 

ql<9cj (^|3y^8c2(n-2)(4nln2n)i_ 

Now taking no larger if necessary, as a > 1, one can ensure that for n > uq, 
{An)i < n, (ln2n)^ < ln(2~^nln2n), 

so 

8c2(n - 2)(4nhi2n)^ < 8c2(n - 2)nln(2"^nln2n). 

Therefore 

q2 ^ g^2 ^^Vg8c2n(n-2)ln(2-lnln2n) 



< l^^y (^6^(^-2) ln(2-lnln2n))^ 



4c2n 



Finahy by ()12p we obtain 

< 94 



, 2 / \ 4c2n 



< C 



with C = 9cf , c = C3 and 7 = 2 + 4c2 . This ends the proof. □ 



4 Proof of Theorem 12.3 

The proof of Theorem 12.31 will be presented in section 14.21 but first in sec- 
tion 14.11 following |LM05] , we will explain how the mechanism of instability 
that we explained in the Gevrey context can be (partly) generalized to an 
analytic context. 

Recall that the first feature of the mechanism is to study perturbations 
of integrable maps and to obtain a result for perturbations of integrable 
flows by a "quantitative" suspension argument. In the Gevrey case, this 
was particularly easy using compactly-supported functions. In the analytic 
case, this is more difficult but such a result exists, and here we will use a 
version due to Kuskin and Poschel ( |KP94j ). 

The second and main feature of the mechanism is the use of a coupling 
lemma, which enables us to embed a low-dimensional map having unstable 
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orbits into a multi-dimensional near-integrable map. In the Gevrey case, 
we simply used the family of maps -0^ : A ^ A defined as in ([1]) and the 
difficult part was the choice of the coupling, where we made an important 
use of the existence of compactly-supported functions. We do not know if 
this approach can be easily extended to the analytic case. However, by a 
result of Lochak and Marco ( |LM05| ) , one can still follow this path by using 
instead a suitable family of maps J^g : A'^ ^ having a well-controlled 
unstable orbit. 



4.1 The modified mechanism 



1. So let us describe this family of maps : A^ A^, g E N*. We fix 
a width of analyticity o" > (to be chosen small enough in Proposition 14.11 
below). For q large enough, Tg will appear as a perturbation of the following 
a priori unstable map 



J7 = ^f(/2+/2)+cos27r9i . ^2 



A^ 



More precisely, for q € N*, let us define an analytic function /g : A^ — )• M, 
depending only on the angle variables, by 

where 

= (sin7rei)-(«'-), = -7r-i(2 + sin27r(02 + 6-i)). 

We still have to define the exponent z^(g, a) in the above expression for fg^\ 
Let us denote by [ . ] the integer part of a real number. Given cr > 0, let Qo- 
be the smallest positive integer such that 



In 



then we set 



v{q,a) 



Ana 
Ing, 



+ 1 



1, 



Aira 



+ 1 



> 



In particular, for q > q^, v{q, a) > 2 and it is always an even integer, hence 
/(^^ is a well-defined 1-periodic function. The reasons for the choice of this 
function f^'^^ are explained at length in |Mar05j and |LM05j . so we refer to 
these papers for some motivations. 
Finally, for g > ^o-, we can define 



jr = o : A2 ^ A^. 



(13) 
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Let us also define the family of points {Cq,k)kez of by their coordinates 

e^.fc : (^1 = 1/2, /i = 2,02 = 0,l2 = q'\k + 1)). 

Clearly, the /2-component of the point ^q^k converges to ±00 when k goes 
to ±00, hence the sequence {Cg,k)kez is wandering in A^. 

2. The following result was proved in |LM05] . Proposition 2.1. 

Proposition 4.1 (Lochak-Marco). There exist a width a > 0, an integer 
qo and a constant < d < 1 such that for any q > qo, the diffeomorphism 
J-q'. A? ^ A? has a point C,q £ A^ which satisfies 

\:F^,'{Q-Uk\<d''^'''''^- 

As a consequence, the orbit of the point Cg ^ A^ under the map Fq is 
also wandering in A^. In particular, for A; = and k = 3q the above estimate 
yields 

\Q - e,,ol < l-Ff '(Cq) - ^,,3,1 < d''^'''^\ 

and as 

\Cq,3q - Cq,o\ = 3, 

one obtains 

|-Ff'(C,)-C,l>3-2d'^(«''^)>l. (14) 

The proof of the above proposition is rather difficult and it would be 
too long to explain it. We just mention that crucial ingredients are on the 
one hand a conjugacy result for normally hyperbolic manifolds (in the spirit 
of Sternberg) adapted to this analytic and symplectic context, and on the 
other hand the classical method of correctly aligned windows introduced by 
Easton. 

3. Now this family of maps : A^ — > A^ will be used in the coupling 
lemma. More precisely, recalling the notations of the coupling lemma 13. H 
in the following we shall take m = 2, 

which is just a rescaled version of the map F^ we introduced before, and 
f = fn = Qn^fqni for some positive integer parameters Nn and g„ to be 
defined below. 

It remains to choose the dynamics on the second factor, and here it will 
be an easy task. In order to have a result for a continuous system with n 
degrees of freedom, we set m' = n — 3, and it will be just fine to take 
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If {pj)j>o is the ordered sequence of prime numbers, now we let Nn be the 
product of the n — 3 prime numbers {pn+A, ■ ■ ■ ,P2n}, that is 

Nn = Pn+iPn+5 ■ ■ ■ P2n- (15) 

The next proposition is the analytic analogue of Proposition \'6.2\ and its 
proof is even simpler. 

Proposition 4.2. Let n > 4 and a > 0. Then there exist a function 
Qn € Aai^'^~^) and a point a„ G A"~^ such an is Nn-periodic for Gn and 
{gn,Gn,an, Nn) Satisfy the synchronization conditions 

gn{an) = 1, dgn{an) = 0, gn{Gn{an)) = 0, dgn{Gn{an)) = 0, 

for 1 < k < Nn — 1. Moreover, there exists a positive constant c depending 
only on a such that if 

g„ = 2Ar„4[e^("-3)P-], (16) 

the estimate 

qn'%n\a<N-^ (17) 

holds true. 

The function gn belongs to >4.(j(T"~^), but it can also be considered as a 
function in Aa{^^~^) depending only on the angle variables. 

As in the previous section, one can easily see that the coupling lemma, 
together with both Proposition 14.11 and Proposition 14.21 already give us a 
result of instability for a perturbation of an integrable map, but we shall 
not state it. 

Proof. Recall that for p G N* , we have defined in Lemma 13.51 an analytic 
function r/p : T — )■ M by 

rjp{9) = -^cos27r/^ . 
\^ z=o / 

We choose our function gn G Aa{T"'^'^) of the form 

gn = gi^^ ^ ■ ■ ■ g^n~^\ 

where 

gliHOi)=r]p„^r+M), 3<i<n-l, 
and our point a„ = {an \ • • • , ai" G A""'^ where 



4*) = (0,p;|i+,), 3<^<n-l. 
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Recalling the definition of G„ and Nn, it is obvious that a„ is A^„-periodic 
for Gn- Moreover, by Lemma |3.5^ the function ry^ satisfies 

r/p(0) = 1, r?;(0) = 0, r/p(fc/p) = r?;(A;/p) = 0, 

for 1 < /c < p — 1, from which one can easily deduce that {gn,Gn,an, Nn) 
satisfy the synchronization conditions ([S]). 
Concerning the estimate, first note that 

\Vp\a < e'^^P 

so that 

\a \ < \n \ ■■■\n I < p47ra(n-3)p2n 

Therefore, if we set c = Sna, then by definition of g„ one has 

> /v2p47ro-(n-3)p2n 
yn — n 

and this eventually gives us 

qn'%n\a<N-^, 

which is the desired estimate. □ 
4.2 Proof of Theorem [2731 

4. First we shall recall the following result of Kuksin-Poschel ( |KP94j . see 
also |Kuk93p . 

Proposition 4.3 (Kuksin-Poschel). Let '■ A"^-*^ A"^-*^ be a bounded 
real-analytic exact- symplectic dijjeomorphism, which has a bounded holo- 
morphic extension to some complex neighbourhood Vq, for some width q> Q 
independent of n £ N*. Assume also that |^'„ — goes to zero when n 
goes to infinity, where h{Ii, . . . , In-i) = + ' ' ' + ^n-i)- 

Then there exist uq G N* , p < g such that for any n > uq, there exists 
fn G ApiT^ X B), independent of the variable In, such that if 

Hn{e, I) = + . . . + Z2_i) + /„ + fn{0, 1), {0, 1) G A", 

for any energy e G M, the Poincare map induced by the Hamiltonian flow 
of Hn on the section {On = 0} H H~^{e) coincides with "^n- Moreover, the 
estimate 

|^'„-^%< |/n|p<5n|^'n-^%, (18) 

holds true for some constant 5n that may depends on n G N* . 
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This suspension result is slightly less accurate (since more difficult) than 
Proposition 13.61 there is a constant 5n depending on n. However, what 
really matters is that the resulting width of analyticity p depends only on g 
and R, but not on n. 

5. Now we can finally prove the theorem. 

Proof of Theorem \2.3[ Let n > 4, R > I and a > given by the Proposi- 
tion and let A^^^ and qn defined as in (fT5]l and (fT6]l respectively. 

We will first construct a map ^„ with a well-controlled wandering point. 
To this end, by Proposition 14.21 we can apply the coupling lemma 13.11 with 
the following data: 

and with the function g„ and the point a„ given by the aforementioned 
proposition. This gives us the following: if 

where ^(^i) = cos27r0i, then the A^^^i-iterates of the map 

= o $'^+^" : A""^ ^ A""^ 
satisfies the following relation: 

^^(x,a„) = ($«"'^^" oF„^"(x),a„), x € A\ (19) 
Now let us look at the map 

^Qn^f.n o F^" = $9"'/9n q 5 +^1 ) + ' COS 27rei ^ ^" ^ 

If Sn{Oi, 02, h, h) = (^ii 021 Nnh, Nnh) is the rescaling by Nn in the action 
components, one sees that 

where J'j^f-^q is defined in ()13p . Now by Proposition 14.11 choosing n large 
enough so that NZ'^Qn > qo, this map has a wandering point Caf-i^ £ A^, 
which by ()14p satisfies 

Using the above conjugacy relation, one finds that the point 

Xn = 5-l(C^-V)GA2 
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wanders under the iteration of •J*'^" o F^^ , and that its drift is bigger 
than one after Nn{3N~'^q'^) = 3N~^q'^ iterations, that is 



1„2 



> 1. 



By the relation ([19]) this gives a wandering point x„ = (xn>«n) S A" ^ for 
the map satisfying the estimate 

|^3'?^(x„)-x„,| > 1. (20) 

Next let us estimate the distance between ^'„, and the integrable dif- 
feomorphism <I>^. First note that since Un,Vn G AaC^^^^), extends 
holomorphically to a complex neighbourhood of size a. Let us now estimate 
the norms of n„, and Vn- Obviously, one has 

iV„~' < \Vn\a < e'^''N-\ 

By definition of fq and the exponent v{q, a), one easily obtains 

U-lf I <o-V2|/'(2)| 
\'in Jqn\cr ^ Hn \J k) 

and hence 

where the last inequality follows from the estimate (jl7p . Then by using 
Cauchy estimates and general inequalities on time-one maps, we obtain 

N-^ < l^n - < c,iV„-', (21) 

for n large enough, and for some constants Ca and g > depending only on 
a (for instance, one can choose g = 6~^a). 

Now we can eventually apply Proposition 14. 3t there exist no € N*, p < 
such that for any n > uq, there exists fn G ^pCT" x B), independent of the 
variable In, such that if 

Hni9,I) = ^ilf + ---+lLl)+In + fn{e,I), (6 , 1) G , 

for any energy e G M, the Poincare map induced by the Hamiltonian flow 
of Hn on the section {On = 0} n H~^{e) coincides with Clearly, the 
wandering point Xn for ^n gives us a wandering orbit {x{t),t, In{t)) = 
{x{t),9n{t), Init)) for the Hamiltonian vector field generated by Hn, such 
that 

x{k) = ^'^{xn), kez. 

In particular, after a time t„ = 3g^, by the above equality and the rela- 
tion (1201) this orbit drifts from to 1. 
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Now it remains to estimate the size of the perturbation = \fn\p and 
the time of drift t„ in terms of the number of degrees of freedom n. First, 
by dlH]) and 

< En < CnN;;\ (22) 

with Cn = C(j5n- Then, by the prime number theorem, taking uq large 
enough, one can ensure that 

^'2n/4 < Pn+i < P2n, 4 < i < n, 

which gives 

(p2n/4)"-3 <Nn< iV„^ < P2n < 4iV„^ . (23) 

We can also assume by the prime number theorem that for n > hq, one has 
2n In 2n < p2n < 2(2n In 2n) = 4n In 2n. (24) 
From the above estimates (|23p and ()24p one easily obtains 

g(ra-3)ln(2-inln2n) ^ ^ ^ g(n-3) ln(4n In 2n) ^25) 

and, together with one finds 

g-2(n-3)ln(4nln2n) ^ ^ ^ ^ g-2(n-3) ln(2-ln In 2n) ^26) 

Concerning the time r„, we have 

1 

Tn = Ml < I2iv8e2^('^-3)P2„ < i2iv8e8'=("-3)^"'^ , 
where the last inequality follows from ()23p . Then using ()25p we have 



1 



A^^"-'' < 4n In 2n 

and from (1221) we know that 



-"n — 



, 4 



Sr. 



SO we obtain ^ 

^2 ^ 19 / "^3 A 32c(n-3)nln2n 

Then one can ensure that for n > no, 

ln2n < ln(2"^nln2n), 

so 

32c(n - 3)nln2n < 32c(n - 3)nln(2"inln2n). 
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Therefore 

g2 <; 12 I — ] g32c(n-3)nln(2-lnln2n) 



, 4 



Sr. 



Finally by (|26p we obtain 



7^ 

in 



4 / \ 16cn 

< 12 I £!i ] ( £!i 



with C = 12 and 7 = 4 + 16c. This concludes the proof. □ 

A Gevrey functions 

In this very short appendix, we recall some facts about Gevrey functions 
that we used in the text. We refer to [MS02], Appendix A, for more details. 

The most important property of a-Gevrey functions is the existence, for 
a > 1, of bump functions. 

Lemma A.l. Let a > 1 and L > 0. There exists a non-negative 1-periodic 
function (pa,L £ G"'^ ^]) whose support is included in [—\, j] and such 
that ipa,L{0) = 1 and ip'^ ^(0) = 0. 

The following estimate on the product of Gevrey functions follows easily 
from the Leibniz formula. 

Lemma A.2. Let L > 0, and f,g e G"''^(T" x B). Then 

\fg\a,L < \f\a,L\g\a,L- 

Finally, estimates on the composition of Gevrey functions are much more 
difficult (see Proposition A.l in |MS02] ). but here we shall only need the 
following statement. 

Lemma A. 3. Let a > 1, Ai > 0,Li > 0, and L,J be compact intervals of 
M. Let f e G"'^i(/), g G G°''^^{J) and assume g{J) C /. // 

\9\aM < 

then f og £ G°'-^i (J) and 

\f og\aM ^ l/U,Ai- 
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